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Abstract 

We study systems containing electrons and nuclei. Based on the fact that the 
Thermodynamic limit exists for systems with Dirichlet boundary conditions, we 
prove that the same limit is obtained if one imposes other boundary conditions such 
as Neumann, periodic, or elastic boundary conditions. The result is proven for all 
limiting sequences of domains which are obtained by scaling a bounded open set, 
with smooth boundary, except for isolated edges and corners. 

1 Introduction 

We consider systems composed of electrons and nuclei, i.e., point particles which inter- 
act via Coulomb interaction with the negatively charged particles being fermions. Due 
to their important role in describing nature, such systems have been intensively investi- 
gated. In particular, the thermodynamic limit, i.e., the limit in which the system becomes 
large, has been studied extensively in pQ. In that work, it was shown that the thermody- 
namic limit exists for thermodynamic quantities, such as the pressure and the free energy 
density, provided that they are denned using Dirichlet boundary conditions. Further- 
more, it was shown that these quantities possess the properties which are expected from 
phenomenological thermodynamics. 

In order to define the canonical and the grand canonical partition function, one has 
to confine the particles of the system to lie in a bounded set AcR 3 , which we choose 
to be open. For the confined system to be well defined its Hamiltonian should be self 
adjoint. This requires that one imposes suitable boundary conditions on the boundary of 
A. For each particular choice of boundary conditions one obtains a canonical and a grand 
canonical partition function. In order to study the thermodynamic limit one considers 
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a sequence {A;} of bounded open domains such that the volume of A; tends to infinity 
as I — > oo. For systems with Dirichlet boundary conditions it was shown in pQ that 
the canonical and grand canonical partition function exist for a large class of limiting 
sequences {A;}. Moreover, the limit is independent of the particular sequence. 

In this work we prove that, indeed, the same limit is obtained for systems with Neu- 
mann, periodic, or reflecting boundary conditions. We prove our result for limiting se- 
quences which are obtained by scaling a bounded open set, which has a smooth boundary, 
except for isolated edges and corners. This class of limiting sequences is smaller than the 
class for which the thermodynamic limit for systems with Dirichlet boundary conditions 
has been shown to exist. We want to point out that this is only partially technical. For in- 
stance, there exist sequences of domains for which the thermodynamic limit of the ground 
state energy for Dirichlet boundary conditions exists, whereas for Neumann boundary con- 
ditions the ground state energy diverges to — oo. Although such sequences are somewhat 
pathological, this demonstrates that the independence of boundary conditions for systems 
composed of electrons and nuclei cannot be considered as trivial. We will also comment on 
possible more general classes of limiting sequences for which our proof is applicable. For 
notational simplicity, we only state and prove our results for systems composed of a single 
species of negatively charged fermions and a single species of positively charged particles 
being bosons. The results as well as their proofs generalize to multicomponent systems 
in a straight forward way. We state the main result and present its proof for both: zero 
temperature and nonnegative temperature. Despite that the latter implies the former, 
we present that way an independent and technically easier proof for the temperature zero 
case. 

To prove the independence of the boundary conditions we use a sliding technique, 
which was introduced in |2j, and refined in Thereby, one decomposes the space 
into simplices. By sliding and rotating the simplices one obtains a lower bound for the 
Hamiltonian of a large system in terms of Hamiltonians defined on the smaller simplices. 
Simplices which lie in the interior of the large system have Dirichlet boundary conditions. 
Whereas simplices on the boundary, i.e., simplices which intersect with the boundary of 
the large system, are subject to mixed boundary conditions. Using that the many body 
Coulomb potential can be estimated below by a sum of one body potentials [8 j , we then 
show that the thermodynamic quantities in the boundary simplices are bounded. In the 
thermodynamic limit the sum of all the boundary contributions is proportional to the 
surface. This is negligible compared to the bulk contribution, which is proportional to 
the volume. 

We want to point out that independence of boundary conditions has been studied for 
systems with hard core interactions (see [3], [5], and references given therein). 

The paper is organized as follows. In Section El we introduce the model and state the 
results. In section |3] we present the proofs. 

2 Model and Statement of Results 

We shall first recall the definition of Dirichlet and Neumann boundary conditions 0. 
Let A be a bounded open set in IR 3 . The Dirichlet Laplacian for A, — A°, is the unique 
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self-adjoint operator on L 2 (A) whose quadratic form is the closure of the form 



\W<p\ 2 dx 

'A 

with domain Cg°(A). The Neumann Laplacian for A, — A^J, is the unique self-adjoint 
operator on L 2 (A) whose quadratic form is 

\V(p\ 2 dx 

J A 



with domain -f^(A) = { / G L 2 (A) | V/ G £ 2 (A) (in sense of distributions) }. 

The model consists of electrons {h 2 /2 — 1, m — 1, |e| = 1) and nuclei with mass M 
and charge z. We assume z to be rational. The electrons are fermions, while the statistics 
of the nuclei is irrelevant. Let A e M 3 be an open set. The Hilbert space Hn,a, with 
N = (n, k) GN 2 , for n electrons and k nuclei is the subspace of L 2 (A xZ 2 f n ® L 2 (A) m 
carrying the permutation symmetry appropriate to the given statistics. The Hamiltonian, 
acting on Hn,a, is 
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where the electron coordinates are x«, the nuclear coordinates are Ri and by B we denote 
the type of boundary conditions, e.g. N, D, and M stands for Neumann, Dirichlet, and 
mixed boundary conditions. Variable particle numbers are accounted for by means of the 
direct sum 

= (J) Kn,a , Hj* = (J) H% A . (1) 

N N 

The grand canonical partition function and the (finite volume) pressure are defined by 

z B (p, n, A) = Tr nA e -^-^ = £ Tr nNA e -/W. A -M-~) 

/(AM,A) = (/^lArMogS^^^A), 

where = (//„, G M 2 stands for the chemical potentials of the electrons and the nuclei, 
(3 > is the inverse temperature, and N denotes the particle number operator, for which 
we use the same symbol as for its eigenvalues. Here and below the volume of a subset 
Q in M 3 is denoted by The canonical partition function of the system at reciprocal 
temperature (3 and the free energy per unit volume are defined by 

Z B (P,N,A) = Tr^ A e-^ 
f B (P,N,A) = -([3\ AD" 1 log Z B (f3,N, A). 
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Furthermore, we consider the following zero temperature expressions, which we will denote 

as 

G B (p, A) = inf a HK {H* — p ■ N) , g B (p, A) = j^G B (p, A) , 

E B (N, A) = inf a nA (H B tA ) , e B (N, A) = ±-^E B (N,A) . 

Definition 1. A sequence {A;} of bounded open sets in IR 3 is called a regular sequence of 
domains if: 

(i) For I — > oo ; |A;| — > oo. 

(mJ For eac/i fixed h > as I —> oo f wtt A^ = IR 3 \ A^ 

|{ xG A ( | d(x,kf) < h }|/|Aj| — > and |{ x G Af | d(x, A,) < /i }|/|A,| -> . 

(Hi) There exists a 5 > swca that for all I, \Ai\/\Bi\ > 5, where Bi is the ball of smallest 
radius containing A/. 



It was shown in |T| that for regular sequences {A/} the thermodynamic limits 
p D ((3, p) = lim p D (/3, p, A,) , lim g D (p, A,) = g D (p) 

l—yoo I— >oo 

for Dirichlet boundary conditions exist and are independent of the particular sequence. 
To study the thermodynamic limit for the canonical ensemble, we consider systems with 
no net charge, i.e., 

N = (n, k) with k = zn . 

We introduce the set 

PS = { (Pe, Pk, ) € I Pk = Zp e } , 

corresponding to neutral charge configurations. In [1 , it was also shown that for a regular 
sequence of domains {A/} and neutral {N{\, i.e., Ni G N 2 fl Ps, with 

Ni 

, lim u~7 = P = <^ e ' ^ fc ) e P s ' 

i— >oc |A|| 

the limits for Dirichlet conditions 

lim f D ((3, AT,, A,) = / D (/?,p) , lim e D (N l ,A l ) = e D (p) 

(— +oo I— >oo 

exist independent of the particular sequence and are convex functions of p G Ps- The 
value of p gives the density of the particles. Furthermore, it was shown that the canonical 
and the grand canonical ensembles are equivalent, i.e., that 

p D (P, p) = sup (p ■ p - A/3, p)) , g D (p) = inf pePs (e D (p) - p ■ p) . (2) 



We note that (J2J) implies that g D (p) is concave and that p D (/3,p) is convex, and hence 
they are continuous functions of p. We now state our main result. 
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Theorem 1. Let A C M 3 be a bounded open set with smooth boundary, except for isolated 
edges and corners. Consider the sequence A L = LA for L > 0. Then 

(a) £ N (p) = limx^oo |A L |- 1 G N (/x ) A L ) = £ D (p). 

(b) p N (/3,p) = \im L ^ ooP N ((3^,A L )=p D ((3, f i). 

We want to point out that sequences satisfying the assumption of Theorem [T] are regular 
sequences of domains. Theorem ^ has the following consequence. 

Corollary 2. Let {A^} be a sequence of domains as in Theorem^ Let {Nl} be a 

sequence with neutral charge configuration, i.e., G N 2 fl Ps, such that 

r hm ta — r = p G Ps ■ 

Then 

(a) e N (p) := lim^ e N (AT L , A L ) = e D (p) ; 

(b) / N (/3, p) := lim^ / N (/3, JV L , A L ) = / D (A p) • 

Remark 1. Note that we only consider systems which consist of electrons and one type 
of spinless nuclei. The results and their proofs generalize in a straight forward way 
to multicomponent systems, with all negatively (or positively) charged particles being 
fermions. 

Remark 2. The essential technical requirement in the proof of Theorem^on the sequence 
of domains A^, apart from being regular, is that the thermodynamic quantities of the 
boundary simplices are bounded, cp. Lemma |SJ This in turn holds for all sequences 
satisfying the assertion of Lemma which is stated in the next section. We want to point 
out that there do exist sequences of domains for which the thermodynamic limit does not 
exist for systems with Neumann boundary conditions, and yet for Dirichlet conditions the 
thermodynamic limit exists, cp. the example below. 

Example 1. We consider a system where the charge z of the nuclei is one. Let {A;} be 
the union of a large ball B\ of radius / and a shrinking ball Bi-a of radius /~ 4 separated 
from the large ball by a constant distance. Let {N{\ be a sequence with Ni G N 2 n Ps 
lim^oo |Ai| _1 A/"; = p. The sequence {A;} is a regular sequence of domains. We place 
one electron and a single nucleus in the small ball and put both in the Neumann ground 
state. In that situation the small ball has neutral charge distribution and hence there is 
no Coulomb interaction with the large ball. This provides us with the following upper 
bound 

e N (A ; , Ni) < \A l \- 1 E D (B l , N t - (1, 1)) + lA^I^H -2 / T—^l dxd V ■ 

i(B i _ 4 ) 2 F - y\ 

The first term on the right hand side converges to e D (p) while the second term diverges 
to — oo. The same conclusion is easily seen to hold if we connect the small ball to the 
large ball by a thin tube provided that its thickness shrinks fast enough. 
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Finally we want to consider more general boundary conditions. Consider for instance, 
a Laplacian — with boundary conditions such that 

-A* < -Ai < -A? . (3) 

(Here and below operator inequalities are understood in the sense of forms [6 .) Then 
Theorem^ and Corollary respectively, imply that the same limits are obtained for sys- 
tems with boundary conditions satisfying (|3*|). We note that periodic boundary conditions 
are of this type. 

Elastic boundary conditions, with elasticity a, are defined as follows. Let t a k denote 
the quadratic form which is the closure of the form 

|V0| 2 dx + a / \<p\ 2 dS 

J A Jdk 

with domain H 1 (A) D C(A) and where dS is the surface measure of dA, the boundary of 
A. Let — A^ denote the unique self adjoint operator with quadratic form t a A . Functions 
in the domain of — A^ satisfy 




at the boundary of A, where d<p/dn denotes the normal derivative. Note that boundary 
conditions with elasticity zero are Neumann boundary conditions. For positive elasticity, 
a > 0, we have the operator inequality — A^J < — A^ < — A^. This implies the statement 
of the following theorem in the case where the elasticity is positive. That it, indeed, holds 
for negative elasticity will be shown in Section 13.41 

Theorem 3. Let Al be as in Theorem^ Then, for elastic boundary conditions with real 
elasticity a, limL_ >00 p CT ( ( 5, fj,, Al) =p D ((3,n). Let {Nl} be a sequence as in Corollary^ 
Then lim^ / CT (P, N L , A L ) = f D ((3, p). 



3 Proofs 

First we show that Corollary El follows from Theorem ^ In subsection 13.21 we will prove 
Theorem Q which is our main result. The prove is based on a lemma which estimates the 
contributions from the boundary terms. The prove of that lemma is deferred to subsection 
13.31 In subsection 13.41 we prove Theorem El concerning reflecting boundary conditions. 

3.1 Proof of Corollary [2] 

(a). We know that for p G P$ 

e D (p) >liminfe N (AT i ,A L ) . (4) 

For a given p £ P$ there exists, by the convexity of e D , a p such that 
e D (p')>e D (p) + p-(p'-p) , for all p' e P s . 
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Hence 

inf p / ePs (e D (p) - p ■ p) = e D (p) - p ■ p . (5) 

We have 

liminfe N (iV L ,A L ) = liminf Ai) 



L^oo L— >oo hVi> 

. t E"{N L ,A L )-n-N L , 
= Iim mi — — : \- u, • p 

L^oo \A L \ 

. r . f /. f £ N (iV,A L )-p.A^ 

> hm mi mi /v n — i + M • P 

L-oo V |A L | J 

= Ap) + P"P 

= inf p , ePs (e D 0') - /x- p') + p- p 
= Ap), 

where we have used Theorem Q (a) in the fourth, eq. (J2J in the fifth and eq. (0) in the 
last line. The above inequality together with proves (a). 

(b). The proof of (b) is analogous to (a). We know that for p E Ps 

-f D ((3,p) <limsup(-/ N (/?,AT L ,A L )) , (6) 

where we used that the map A \— > Tr e A is operator monotone. Since Pg 3 p i— > — / D (/3, p) 
is concave, there exists for a given p £ P$ & p such that 

-f(P,p')<-f((3,p)-p,-(p'-p) , for all p' G P 5 , 

and hence 

sup (p' • p - A/3, P')) = P ■ P - A/?, P) • (7) 
p'eP s 

We have 

limsup (-f N (P,N L ,A L )) = 

L— »oo 

= limsup (-p • p + OSlAil)- 1 /?/* ■ N L - f N ((3, N L , A L )) 
= -p ■ p + limsup ((/^l)- 1 logZ N (/3, N L , A L )e^ N ^) 



L— >oo 



< -p • p + limsup ( (/^lAil)- 1 log ( ^ N (A ^ A i) e/3/ " 

L^oo y y N 



N 



= -p-p + p D (/3,p) 
= -f D (P,p), 

where we have used Theorem ^ (b) in the fourth and eqns. ()2I7|) in the last line. The 
above inequality together with (jBJ) proves (b). □ 
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3.2 Proof of Theorem Q] 

To prove Theorem^ we will make use of the localization method of [3j (see also [2]), where 
one breaks up R 3 into simplices in the following way. Cutting the unit cube W = [0, l] 3 
with all planes passing through the centre and an edge or a face diagonal of W, one 
obtains congruent simplices A n C W, (n = 1, ...,24). The simplices A a = A n + z, with 
a = (z,n) G Z 3 x {1,...,24} =: I, yield a partition of M 3 up to their boundaries. We 
then choose a spherically symmetric <fo G 0^°(IR 3 ) with J ip\ = 1 and {<^ (x) 7^ 0} = 
{\x\ < 1/2}. Let Xa be the characteristic function of A a . Setting (p(x) = r]~ 3 / 2 Lp (x/ri) 
and j a = (Xa * V^ 2 ) 1 ^ 2 ) we obtain a partition of unity, i.e., 



1 , (x G 1R 3 ) 



with j Q G C°°(1R 3 ). There are congruent simplices A+, which are scaled copies of A Q , 
such that 

supp j a C A+ 
|AJ|<|A a |(l + 0(^)) 

as 77 I 0. 

The following definitions depend on rj and I > although the notation will not reflect 
this for simplicity. For the moment, let A C I 3 be any bounded open set. For y G W and 
R G 50(3) we set 

A y,R = R -i A _ i y _ 

We define the set 

1(A) = { a e J|/A+ n A y ' R ^ for some yeW,Re SO(3)} . 

For a G 1(A), let 7i a := TC lA + be the many particle space for the simplices ZA+ as given 
in (^). By H a A yRA we denote the Hamiltonian on ZA+ R A 2 '' with Neumann conditions 
on (ZA+) fl dA y,R and Dirichlet conditions on the remaining part of the boundary. The 
operator H^f y R A acts on 

and hence on TC a via the canonical embedding. Note, if / A+ C A y ' R , then H^ yRA = H? + 
and Ti,a,y,R,A = 7~t a . We define the Hilbert space and a Hamiltonian, acting on it, as the 
direct integrals 

Hi(A) = / dydfi(R) 6?) H a , 

7iy X 50(3) aeJ(A) 

^/ N (A) = / W*) E ^.fl-A > 

7iyx50(3) Q£/(A) 

where c?// denotes the Haar measure on SO (3). We shall define a map J : H\ — > l~ti(A) as 
follows. Let j y ,R t a '■ L 2 (A) — > L 2 (ZA+) be given by 

{jy,R,ai>){x) = j a (x/l)i>(R(x + ly)) . 
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Define 

Jy,R ■■ L 2 (A) 

jy,R 

This lifts to a map between the many particle spaces 

r(j y , R ) :H A ^ (g) n a , 

ae/(A) 

which acts as the iV-fold tensor product of j y ,R on iV-particle states. We may now define 

r® 

J :H A ^ Hi(a) , J = dyd/j,(R)T(jy iR ) . 

JWxSO(3) 

We note that the map 3* y ^ a '■ L 2 (IA+) — > L 2 (A) is given by 

(f y ,R,ai>)(x) = UR-\x/l) - yMR-'x - ly) . 

Hence jy )R j y ,R '■ L 2 (A) — > L 2 (A) acts as multiplication by J2 a ei(A) ^(-^"H^/O ~~ v)- This 
function of x G A equals 1. We conclude that J* J = 1, i.e., that J is an isometry. We 
state the following Lemma, c.f. Lemma 7 in 

Lemma 4. Let t] = Then 

H]* > nJ*Hf {A) J- /^const • N 
for large I, where < k < 1 and k — 1 + 0(/ _1 ) as I — > oo. 

For the proof of this Lemma we refer the reader to the proof of Lemma 7 in [Hj , where 
the statement is for the Dirichlet Laplacian. With little modification of the proof given 
there one can proof Lemma HJ 

From now on, let A be a fixed open set as in the assumption of Theorem ^ i.e., 
bounded with smooth boundary, except for isolated edges and corners. We will consider 
the sequence of scaled copies A^ = LA, with L > 0. Note that the claim of Lemma |U 
of course holds for all A^. Let A denote a simplex which is similar to one (and thus 
all) of the simplices A a , i.e., equal up to dilations, translations, and rotations. By A c 
we denote its intersection with A L , i.e., A c = A n A L . Let -A^ c (H^J denote the 
Laplacian (Hamiltonian) on A c with Neumann conditions on A D OAl and Dirichlet 
conditions on the rest of <9A C . We note that — A^ c is the unique self adjoint operator on 
L 2 (A C ) whose quadratic form is the closure of the form i— > j A \V<p\ 2 dx with domain 

{ G Lf 1 (A c ) fl C(A C ) | 4> vanishes in a neighborhood of dA n A L }. 

The following Lemma whose prove will be postponed to subsection 13.31 provides us 
with a bound for the contributions from the boundary simplices. 

Lemma 5. Let A^ be a sequence of domains as in Theorem^ and let v > 0. Then there 
exists a number L and constants Ce{h,v) and Cs(/3, fi,v) such that 



a£l{A) 

3y,R,a ■ 

ael(A) 
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(a) G u ( f i,A c )>C E (fi,v)>-oo; 

(b) Tr WAc e~^-^ < C s ((3, ^v) < oo , 
for all A c = An A L , with L > L and | A| < v. 

For the proof of Theorem ^ we will also use the following lemma. 

Lemma 6. Let {/j^} be a sequence in R 2 with lim^oo = and let {A/} be a regular 
sequence of domains. Then 

(a) lim^oo \Ai\~ 1 G D (ii l , A/) = liirn^ |A,|- 1 G D (/x, A ; ) = g D (^) ; 

(b) lim^oop ^,^^;) = linw^O^A,) =p D (/3,M) • 

Proof. We use the notation /x z = (// n ;,/%) and e = (e, e). For e > 0, there exists an Z 
such that for all I > / 

£in - £ < /Ltn; < /in + e 

Mfc - e < /ifc/ < /ifc + e • 

For (a), we note that 

[A i [- 1 C I3 (M + e, AO < lA^'G ^,^) < lA^ 1 ^ - e, A,) , V I > l . 

Hence 

g D (» + e) ^liminflAr 1 ^ ^,^) < limsup lA^ 1 ^/^ A,) < <7 D (m - e) , 

and, by the continuity of /x i— > g D (£t), (a) follows. 

For (b), we first note that, by equation @, ^ i— > p°(f3, n) is convex and hence contin- 
uous. In analogy to (a) we have, using that A \— > Tre" 4 is (operator) monotone, 

p D (/3, /x - e, A,) < p D (/3, M/ , A,) < p B ((3, M + e, A,) , V Z > l . 

Hence 

p D (/3,/x-e) < liminfp^/^A,) < limsupp D (/3, A,) <p D (/3,/x + e) , 

and (b) follows by the continuity of \i \— > p D (/3, /x). □ 

Proof of Theorem^ (a) Since — fi ■ N > if]^ — ^ • AT, the inequality 

g B (^)>\imsnpg N (^A L ) (8) 

L^oo 

is obvious. We shall show the inequality liminf^oo g N (/x, Al) > g B (/-i). We introduce 

2V /(A) = / dj/d//(i2) V AT a , 

J W xSO(3) a6J(A) 
10 



where N a denotes the number operator of TC a - Note that J*Njr\\J = N. By LemmaEJ 

H^-ti-N>K.r (Hf (A) - ii x ■ N m ) J , 
where we have set jl l = (l/ft)(/x + / _1 const). We define 

^(A) = {a|JA+cA**} 
#ji(A) = {a|lA+ n 0A*V 0} ■ 
Let \l/ G T~La l be normalized to one and smooth. We observe that 
(*, (H% L - (J, ■ N)*) 

> K f dydn(R)\T(j y , R )* t Yl ftAA - Ar N a )T(j y , R )* 

>k[ dyd»(R)l Yl G D (»iJK)+ E G M (/i,,/A+nA^ 

JWxSO(3) \ s _,u > 



Hence 



liminf |A L |- 1 G N (/x,A i ) > 

L— >oo 



liminf /c / cfod/i(«) f V ■ |ZA+|" 1 G D (A^A+ 



+ E |^cfc(£i,|JA+|)] 



> ft(l + 0(Z -1 )) • \lA + \- l G D (fi h IA 



where we used that all simplices A+ are congruent to a single one, which we denote by 
A + , and in the last inequality we used Lemma El (a) and that both limits 

; im E ¥T = (9) 

-L— M30 * — ' AfJ 

lim V — — = (10) 

are uniform in i? G S"0(3), y G W 7 . We omit a proof of this simple facts. By Lemma El 
(a), the subsequent limit / — > oo yields 

liminf^Ai)^^^). (11) 

L— >oo 

The two inequalities (JEJ) and (fTTjl show the claim. 
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(b) The inequality 

p D (A^)<limsupp N (/3,/x,AL) (12) 

L^oo 

follows from — A N < — A D . The opposite inequality 

limsupp N (A fi,A L ) <p D (/?,M) 

L— *oo 

is seen as follows. We set to = + / 1 const. Let {<fi}iei be an eigenbasis of J*{kHi^\ l ) — 
l*-iNi(k L ))J ■ Then, using LemmalU we have 

— e ~/3{J'Pi,(^H I ( AL )-iJ.i-N I ^ AL - ) )Jipi) 



iei 

< 



(j<Pi, e-^ KHl ^)-^N I{ A L) ) J(p \ 



< Tr-u e -P(K H i(A L )-v>vNi(A L )) 

< [ dydfi(R) TT Tr^e-PWS-iH-N* 

x TT Tr w RA e -^K^, L -^) 



where, in the third line, we used Jensen's inequality with the spectral measure of {kHi<\ l \ — 
to ■ Ni(a l )) for Jifi. Since < k < 1, we have kH^ > n 2 T® + kV = iJ D _ 1 +, where 
denotes the kinetic Energy, V is the Coulomb potential, and the unitary equivalence 
comes from scaling. Note that all the simplices A+ are congruent to a single one A + . By 
Lemma 0(b), we have 

P N (A/x,Al) 
= (/3|Al|) _1 logS N (/?, fi, A L ) 

< (/SIAlI)^ 1 log (s D (/?, ft, K-HA+y^y^ 1 '^^ • C h (k/3, K"Vi, |ZA+|) sup ^ (|/ ^ (Ai)l) 

< sup(|/^(A L )|)|A L |- 1 |/ t - 1 /A + | .p D (/3,^,K-W) 

+ sup(|/yA i )|)(/3|A L |)- 1 -logC s (^, K -V/,|/A + |) ; 

note that S D > 1 and C s > 1. Thus 

limsupp N (/3,M,A L ) < K- 3 (l + 0(r 1 ))p D (/3,^,K- 1 ZA+) , 



where we have used equations ITUj). Using Lemma E](b), the subsequent limit I — > oo 
gives 

hmsupp N (/?,M)<P D (AM). (13) 

L^OO 

The claim in (b) follows from eqns. (jl2j) and (|13p. □ 
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3.3 Proof of Lemma E] 



To prove Lemma El we will first state a technical Lemma reflecting the geometry of A. 
We recall that A denotes a simplex which is similar to one of the A a , i.e., equivalent up 
to dilations, translations and rotations. 

Lemma 7. Let A^ be a sequence of domains as in Theorem^ and let v > 0. Then there 
exist a constant C\ > and a number L ( depending on v) such that for all L > L and 
all simplices A with | A| < v which intersect dA^, we can choose an open set V containing 
A c = A n A L and smooth coordinates 



with the properties: 

(i) B is a ball centered around the origin and AnV corresponds to either the half space 
restriction {x G V \ y${x) < 0} ; the quarter space restriction {x G V \ y 2 {x),y^{x) < 
0}, or the octant restriction {x G V \ yi(x) , y 2 (x) , y 3 (x) < 0}. 

(ii) The Jacobian Dip has determinant one and Dip~ 1 {Dip~ 1 ) T > C\. 

We want to point out that in the case where A is a box, this Lemma follows trivially 
by choosing C\ = 1 and the coordinate maps to be an appropriate composition of a 
translation followed by a rotation. In that case, also the next lemma is trivial. 

Proof. By assumption, A is a bounded subset of R 3 with smooth boundary, except for 
isolated edges or corners. This means that around any point xo G OA on the boundary of A 
there is an open neighborhood V XQ on which we may choose smooth coordinates (yi, y 2 , 2/3) 
such that A PI V Xo corresponds to either the half space restriction {x G V xo \ y 3 (x) < }, 
the quarter space restriction {x G V XQ \ y 2 [x) , y-$ [x] < }, or the octant restriction 
{x G V Xo I ?/i (x), y 2 (x), 2/3 (x) < }. By rescaling a coordinate we can achieve that the 
coordinate map f has Jacobian determinant equal to one 1 . By possibly adjusting the 
coordinates and choosing the neighborhood V XQ smaller, we can achieve that the images 
of the coordinate neighborhoods V xo under the coordinate maps are balls centered at the 
origin. By compactness there exist constants C\ > and r > such that for each point 
Xq G dA we can choose a coordinate map if : V xo — > B such that 



and moreover \x — Xo| < r implies x G V XQ . Given such a collection of coordinate maps 
for A we obtain, by scaling, a collection of coordinate charts for with properties (i) 

lr This can be achieved as follows. Let (f> : x 1— > w(x) = (wi(x), W2{x), 103(2;)) be a coordinate map with 
Jacobian determinant not necessarily equal to one. Then, we define new coordinates 

pW 3 



It follows that dy = \ det D<fi 1 1 dw and thus the Jacobian determinant of the coordinate map x 1— > y(w(x)) 



if : V — > B jh 



vfc) = (yi( x ),y2{x),y 3 {x)) 



Dip-'iDip-'f > C A 




is 1. 
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and (ii). Moreover, the constant r becomes Lr under this scaling. Thus for large L, 
A c = An Al is contained in some coordinate chart. □ 



Let — AH A j denote the Laplacian on p>(A c ) with mixed boundary conditions, i.e., 
<p maps Dirichlet (Neumann) boundaries of A c to Dirichlet (Neumann) boundaries of 

<p(A c ). 

Lemma 8. Let ip be a coordinate map as in Lemma^ Then the map U : L 2 (V) — ► L 2 (B ) 
f i — > f o ip~ Y is unitary. Moreover, on the form domain of — A A 

-A^ c > U*C A (-A™ Ac) )U . 

Proof. Since the Jacobian determinant of ip is one, C/ is unitary. By abuse of notation 
we write f\y) for (Uf)(y) — f o ip>~ l [y). We set g 1 ^ = (Dip^i^Dip^ 1 ) 7 )^. For functions 
/ in the form domain of — A A we write (/, — /) in terms of the y coordinates and 
estimate 

= C A (Uf,-A™ Ao) Uf). 

□ 

Lemma 9. (Lieb-Thirring estimate) Let v > be fixed. There exists a number L 
and a constant Cm (depending on A), such that for all L > L and \A\ < v we have 

|Tr L2(Ac) (-A^ c + V)_| <C M [ \V4x)\ 5 / 2 dx, 

Ja c 

where V is any locally integrable function on A c = A fl with negative part V- E L 5 ^ 2 . 
(Note that Tr A_ denotes the trace over the negative eigenvalues of the self adjoint operator 
A.) 



Proof. We first observe that if A fl = then the estimate is a simple consequence 
of the classical Lieb-Thirring inequality |7] since we have Dirichlet boundary conditions 
on the whole boundary. Thus for a given v let L Q be as in Lemma [7| Assume that A 
intersects with the boundary of A^ and that |A| < v. Let tp : V — ■> B Q be a coordinate 
map with the properties as stated in Lemma Thus A c G V. We consider first the case 
where V H A L = {x E V \ y?,(x) < 0}. On B we define the reflection r : y 2 , y 3 ) i— > 
(yi,y2, —Hz)- By ip(A c ) T we denote the interior of the closure of <p(A c ) Ut(</?(A c )). Given 
a function h on <^(A C ) we extend it to a function h T defined a.e. on ip{A c ) T by setting 



h T {y) 



h(y) , if 2/3 < 
h(r{y)), if 2/3 > - 



14 



This establishes the isometric injection 



j:L 2 (^(A c )) - L 2 (^(A c n 
Thus for any locally integrable function IV we have 

j*(-A° ( A c r + ^ = - A ?(A c ) + w / ? 

where — A^ A w denotes the Dirichlet Laplacian on <£>(A C ) T w.r.t. the Euclidean metric 8 % K 
By the Neumann condition, j maps the domain of — A^ A ^ into the domain of — A^ A w. 
We then conclude using Lemma |H1 

Tr L2(Ac) (-A^ c + V)_ > Tr x2(Ac) (-A^ c +y_)_ 

> C A Tr L2(v{Ac)) (-A^ (Ac) +^V_)_ 

> C A Tr L2(v , (Ac) r)(-A° (Ac)T + C A 1 VT)_ 

> -C A 3/2 C LT [ \Vl{y)f 2 dy 

= -2C A 3/2 C LT / \V.(x)f 2 dx, 

where we made abuse of notation by denoting V_ o by V_. In the step before last 
we used the classical Lieb-Thirring estimate with constant Clt- If A fl V has an edge 
or a corner the proof is essentially the same we just have to perform several reflections, 
which affects the value of the constant in the inequality by at most a factor 8, since in 
that case, we have to consider the volume obtained by reflecting (f(A c ) on all Neumann 
planes. Likewise we have to extend functions defined on (p(A c ). The details are left to 

— 3/2 

the reader. It follows that the Lemma holds for Cm = 16 C A Clt- d 

Proof of Lemma^ Let A be a simplex with |A| < v. Let L be sufficiently large such 
that the assertions of Lemmas [7| and El hold. Consider now A c = A fl for L > L . The 
Coulomb interaction is 

^ n k ^ 

V c (x 1 ,...,x n ,R 1 ,...,R k ) = i r ._ T .i ~ Z J2J2 | T ._ pi 

l<i<j<n 1 «=1 j=l 11 

l<i<j<k 1 Jl 

We introduce the nearest neighbor, or Voronoi, cells {Tj}^ =l defined by 

Tj = {x \ \x — Rj\ < \x — Ri\ for all / ^ j } . 
Furthermore, define the distance Dj of Rj to the boundary of Tj, i.e., 

£)j = dist(i? i , dTj) = - minilRi - Rj\, j ^ 1} , 
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By Theorem 6 in [S], we have the following inequality 

n _ k 

V c {x l) ... 1 x ni R ll ... 1 R k ) > -J2 W ^ + 8 z2 J2 D 7 1 ' ^ 14 ) 

i=l j=l 

where for x in the cell Tj 

W(x) 2Z + 1 



x — Rj\ 



We note that, in the situation considered here, the coordinates Xi and Rj all lie in A c . 
Using inequality (|Hjl. we find 

n 1 ^ 

^,a c _ - fx k k > y^hj - fx k k + -z 2 DJ 1 , 

t=i i=i 

with hi = — A^ c x . — IV(xj) — /i n . The fermion ground state energy of Y^7=i ^ * s bounded 
below by 2 ^\ ej, where e 3 - are the negative eigenvalues of hi. Hence by Lemma|ni(/+(x) = 

max(/(ac),0); 

^N,A C ~ ^n n ~~ fJ'kk 

k 

> -2C M [ \(W + ^ n ) + \^ 2 dx-ii k k + lz 2 y / Dj 1 

Ja c 8 , =l 

> -2"/ 2 C M [ (\wf 2 + \^ n+ f 2 )dx-^ k k + \z 2 J^Dj\ (15) 

Ja c 8 j=1 

where, for the second inequality, we have used that (a+6) 5//2 < 2 3//2 (a 5 / 2 + 6 5 / 2 ) for a, b > 0. 
We estimate the first term using 



k 

j \wf 2 dx = f Wj(x) 5/2 d 

V / (2^ + l) 5 / 2 |a;-i? J |- 5/2 ^ 

+ Y [ (2z + lf' 2 R-V 2 d 



X 

— JT^Ac 
k 

< V I {2z -i 1 )""'1.r - R,\~'"-d.r 

1= iJ\x~Rj\<R 

k 

J=l ^ lx-«jl>ij 

< (2z + l) 5/2 (4vrLR 1/2 + | A c |i?- 5/2 ) 

< (2^+l) 5 / 2 6^y /6 |A| 1/6 fc 5/6 , 

where we have made the optimal choice for R. To estimate the term involving the Dj, we 
note that for k > 2 , 

k 

J>J<A|A|, (16) 
i=i 
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for some constant A > 0. Using Holders inequality, i.e., 

k / k \ 3 / 4 / k \ V4 

3=1 \j=l / \j=l 

we find 

k 

fc 4 /3 A -V3| A |-l/3<^^-l 

for k > 2. Inserting this into (jl5j). we have 

tf^,A e - /^n - ^ > -Ci| Al 1 / 6 ^ 6 - C 2 | A||/i n+ | 5 / 2 - ^ + CaiAr 1 /^ 4 /^! - S kl ) 

for some positive constants < Cj < oo , (z = 1,2,3), which depend only on z, and A. 
The case k = 1 is accounted for by (1 — Ski)- We minimize with respect to k with the 
result that 

h n,Ac ~ Vnn - ^kk > C e (ij,, v) , V I A| < D , 
for some constant Ce{h,v) G K. Hence we have shown (a). 

To show (b) we decompose the kinetic energy 

n k 

r£ = -5>l, ai -(i/M)5>l A 

i=l j=l 

and use the same calculations as in (a). As a result 

H-N,A C ~ l^n n ~ l^kk 

= ^ + 1(^ + 2^ + 2/^ + 2^*;) 
> ^ + 0(|A|,/i,A;), 

with 

0(| A|, M , k) := -2 3 / 2 d| Al 1 / 6 ^ 6 + C 2 | A||2 Mn+ | 5 / 2 - 2/x fc * + C 3 | - 5 kl ) . 
We estimate the grand canonical partition function as follows 
~ M (/3,/x,A c ) = Ti HAc e-^l-^ 

oo 



n=0 

oo 

1 v^fc AM 



x Tr^(A c )®* Ej=1 Aa -^- • e - mA ^ k) . (17) 

fc=0 



If A does not intersect A^ then we have only Dirichlet boundary conditions and in this 
case it is known that the desired bound exists. It remains to consider the case where A 
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intersects with the boundary of A L . Let ip : V — > B be a map with the properties as 
given in Lemma We shall first consider the case where V (1 = {x ^ V \ y^(x) < } . 
We now use the reflection argument and the notation as introduced in the proof of Lemma 
IH1 There we have shown that on the form domain of —A 



M 



-A^ c > U*C A (-A™ Ac) )U = U*j*C A (-A» (Ac)T )jU . 
Using this estimate we find 

oo oo 
n=0 n=0 

The right hand side of this equation is the grand canonical partition function of an ideal 
Fermi gas with Dirichlet boundary conditions, which is known to be bounded above. 
Similarly we estimate 



I M ^ 3/2 



where the last inequality follows from a standard estimate [§]. Note that |<^(A C ) T | < 2| A|. 
We insert the above inequalities into eq. ()17j) . The sum over k converges, due to the term 
with fc 4//3 . Thus we have shown (b) for the case where V fl does not have any edges 
or corners. If V fl Ax has an edge or a corner the proof is essentially the same we just 
have to perform several reflections. We leave the details to the reader. It turns out that 
in the estimates above (p(A c ) T is replaced by the volume obtained when reflecting (f(A c ) 
on all Neumann planes. Each of the three cases gives us a constant. Taking the largest 
we obtain the desired bound. □ 



3.4 Proof of Theorem [SI 

As mentioned in Section |21 the case a > is trivial. Thus let a < 0. Everything in the 
proof of Theorem holds if we replace Neumann boundary conditions with elastic bound- 
ary conditions. The only part of the proof which does not generalize trivially to elastic 
boundary conditions is the proof of Lemma El We circumvent this by showing that the 
Laplacian with elastic boundary conditions can be estimated below in terms of the Lapla- 
cian with Neumann boundary conditions. We recall that — A^ c is the unique self adjoint 
operator on L 2 (A C ) whose quadratic form is the closure of the form i— > J A \V<p\ 2 dx 

with domain V = { (f> e if 1 (A c ) nC(A c ) | <p vanishes in a neighborhood of dA fl A^ }. Let 
— A^ ,<T be the unique self adjoint operator on L 2 (A C ) whose quadratic form is the closure 
of the form 

0h+ f \X7<f)\ 2 dx + a [ \<f)\ 2 dS 
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with domain D. Below we will show that for all A c = An Al, with |A| < v and L > 1, 



-A£T > r(-A&) - C , (« 



for some r, with < r < 1, and some finite constant C > depending only on a and the 
geometry of A. Thus setting c = (C, C) we have 



H A C = + A C + V 

> tT^ c + V-c-N 

-1/_2thM 



> t (r T^ c + tV) — c - N 

= r-'H^-c-N . 

By 

it is now evident that for a < an analog of Lemma El for elastic boundary conditions 
holds. 

It remains to show (JTEJ). Let 

£ : A -> M 3 

X h-> = (£i(x),£ 2 (^),^3(^)) 

be a real vector field continuously differentiable in the closed region A and satisfying 
the boundary condition n ■ £ < —1 on dA where n denotes the inward normal. First 
observe that such a vector field exists. If A is a box or has smooth boundary this is 
clear. Consider now the general case, where the boundary of A has isolated edges and 
corners. Since A is bounded we can cover it with finitely many sufficiently small open sets 
with the property that on each of these sets we can choose coordinates x ^ y(x) = 
(yi(x) , y 2 (x) , y 3 (x)) such that A n corresponds to either V 7 , the half space restriction 
{x e \ y3^x) < 0}, the quarter space restriction {x G V 1 \ y 2 (x),y 3 (x) < 0}, or the 
octant restriction {x G \ yi(x) , y 2 (x) , ys(x) < 0}, and such that there exists a vector 
field on V 1 which is constant in the coordinate chart and satisfies the required property 
on V^. Pasting these local vector fields together by means of a partition of unity on A 
subordinate to the open covering, we obtain a smooth vector field such that n • £ < — 1. 
Given such a vector field on A, then £l(x) = £(x/L), for L G M + , is a vector field on A^ 
with n • £ < —1 (here n denotes the inward normal of A^). For <fi G H 1 (A C ) fl C(A C ) 
vanishing in a neighborhood of dA fl A^ we have 

/ WdS<( \^ L (-ndS)= [ V(£ L |0| 2 )rfx, 

J And A L JAndA L J A c 

where the equality follows from Gauss' Theorem. We calculate 

V(£l|0| 2 ) = (V£ L )|0| 2 + a(V0*)0 + £ L 0*(V0) , 
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and for any e > we have 

|£l(V0*)0 + £l0*(V0)| < I| V 0| 2 + e|a0| 2 . 

e 

As a result 




This implies (jl8j) . □ 
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